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Abstract: Negative dimensional integration method (NDIM) is a technique to deal 
with D-dimensional Feynman loop integrals. Since most of the physical quantities in 
perturbative Quantum Field Theory (pQFT) require the ability of solving them, the 
quicker and easier the method to evaluate them the better. The NDIM is a novel and 
promising technique, ipso facto requiring that we put it to test in different contexts and 
situations and compare the results it yields with those that we already know by other 
well-established methods. It is in this perspective that we consider here the calculation of 
an on-shell two-loop three point function in a massless theory. Surprisingly this approach 
provides twelve non-trivial results in terms of double power series. More astonishing than 
this is the fact that we can show these twelve solutions to be different representations for 
the same well-known single result obtained via other methods. It really comes to us as 
a surprise that the solution for the particular integral we are dealing with is twelvefold 
degenerate. 



Massless Feynman integrals. Negative dimensional integration method. 



Keywords: 



Three-point function 



Contents 



[l| Introduction. 



i 



I On-Shell Two-loop Vertex. 



1^ Conclusion. 







A| The Remaining Solutions. 







1 Introduction. 

Negative dimensional integration mettioci (NDIM) was first devised by Halliday a/ [|T], Q 
witli tlie aim of dealing with Feynman loop integrals in quantum field theory (QFT). It 
combines the powerful concepts of analytic continuation, translational invariance and 
dimensional regularization|P, |, |^ (DREG) in such a way that the intricate positive 
dimensional integration is transformed into negative dimensional integration of polyno- 
mial type. In practical terms, what one needs to do is to solve systems of linear algebraic 
equations and gaussian-like integrals. 

Before plunging deeply in the mire of it, let us briefly make some comments. First 
of all, what do we mean by negative dimensions'! Obviously, at their face value they 
can only be fictional. Nevertheless, in a stretch of our imagination, if we allow ourselves 
just the possibility of their existence, maybe we can bring into fruition something ever 
undreamed of before. So, to begin with, let us mirror our reasoning with that behind 
DREG. There, the all important concept is the parameter D, the space-time dimension, 
allowed to assume complex values. However, this in no way means that we want or will 
even define operations like scalar product |^ in some general D to have reality. The 
reason is quite simple: Physics — and even the common sense — tells us that D is a 
positive integer number. We keep our eyes in this fact. One may think that this way of 
expressing the result of a given Feynman integral in an arbitrary space-time dimension is 
very elegant indeed, but at the end of the day, one always has to look into the real physical 
world, that is, that of positive integer D. Actually, in the whole process of DREG it is 
important to keep D arbitrary because the analyticity properties depend on the space- 
time dimension and, of course, any singularities do so depend on it too. Wilson^] noted, 
very early, that the so-called integration in D-dimensions is not in fact real literally. It 
only seems to be. It only behaves like it is. 

So, we too define an integral in D-dimensions by means of an analytic function and 
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then work with general D, negative values included. We will not attempt to make any 
sense of it or speculate about the D < "world". What we do is just allow for it without 
trying to see how real — or how fictional, for that matter — it is, before doing the 
analytic continuation back to D > 0. In other words, we are not concerned with seeking 
any meaning for a negative dimensional world nor seeking any new Physics. 

Let us sketch the methodology proper. The idea is quite simple: we analytically con- 
tinue the Feynman integral we want to evaluate into D < and solve it there. Then, 
the result we get there we bring it back into the realm of D > 0, by another analytic 
continuation. As we shall see, the whole procedure is much easier to do compared to 
positive dimensional techniques. Up to now, for all the Feynman integrals we have calcu- 
lated using NDIM, the results agree with the ones calculated in the positive dimensional 
regime]^, |^, |^. This includes even some light-cone gauge loop integrals [|TU|, which know- 
ingly are harder to solve with other approaches. 

The outline of our paper is as follows: in Sec. 2 we solve explicitly a two-loop Feynman 
integral entering the two-loop radiative correction to the massless triangle diagram using 
the negative-D approach. Then, we show how to analytically continue the result to 
positive D and give the result for two on-shell external legs. In Sec. 3 we conclude this 
work commenting on some new results we have for massive one-loop and off-shell two-loop 
Feynman integrals. Finally, in the Appendix we list all the other remaining degenerate 
solutions, for the sake of completeness. 



2 On-Shell Two-loop Vertex. 

As we have mentioned earlier, our work here is done within the perpective of checking 
NDIM methodology for D-dimensional Feynman integrals. Therefore, we have chosen as 
our "lab test" for it, the evaluation of an integral pertinent to the Feynman diagram of 
Fig.l. It is a two-loop graph with four particles in the intermediate states, i.e., containing 
four propagators. For simplicity we take the massless case and in order to compare our 
result with the one already known in the literature, we take a particular limit of two 
on-shell external legs. 

Since ours is a choice example, we cannot become too excited about NDIM and its 
seemingly simplicity, yet one can convince himself that the task of solving Feynman loop 
integrals is quite easy in this approach (at least in principle). Of course, we can envisage 
and even anticipate some technical difficulties in other contexts, which is inherent to the 
method, such as the necessity of dealing with multi-indexed summations. 

NDIM is implemented with a few simple steps: Firstly, we calculate gaussian or 
gaussian-type integrals^, which are not difficult to handle, and some books on QFT even 
list them in tables. Then, one makes an expansion in Taylor series of the result obtained, 
and compare it with the expansion in Taylor series of the original gaussian or gaussian- 
type integral. Comparison term by term of both series then yield a system of linear 

^in D-dimensional momentum space 
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algebraic equations, with constraints arising from intermediate multinomial expansions. 
NDIM requires that we solve this system. The principle therefore is quite simple! Let us 
then take a practical example and see how it works. 



P 




t = k — p q — r — k k 



Figure 1: Two-loop three point vertex. 



Let our launching-pad gaussian-like double integral be 



I — j j dPr d^q exp —aq^ — (5{q — pf — ^r^ — uj{q — r — ky 



(2.1) 



This clearly is a pertinent integral to the diagram of Fig.l. For reasons of simplicity 
and future comparison, let us consider that two of the external particles are real, i.e., let 
them be on-shell, namely, k"^ — — 0. 

Completing the square in the variable q we can carry out the first integration and get. 
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2 I /S^pf., 2/3w£fc f n / 9 \ [1/9 

e + X + X J d r exp y-6r - 2urkj exp - \ uj r 



+2l3ujpr + 2u;^rk^ 



(2.2) 



where 6 = j + uj and A = a + P + u. Following the same procedure we perform the 
remaining integration, remembering that t — k — p and that because of the on-shell 
condition — k"^ — 0, p"^ — 2pk: 



TT 



2\ 2 



exp 



-{aP^ + al3uj)p^ 



(2.3) 



where = 0:7 + au + /?7 + (3uj + 7tc;. Expanding the exponential in Taylor series and 
using the multinomial expansion in 0, we get 
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{n~ro} ni!n2!n3!ra4!n5!n6!n7! 

^Q,ni+n2+"-3+"-4^ni+n2+n5+"6^"-i+"3+n5+n7|^"-2+"-4+".6+"-7 ^2 4) 

where the summation indices must satisfy the constraint —rii — 712 — \D = 713 + 77.4 + 71,5 + 
Tig + 777 which is the condition imposed by the multinomial expansion. 

Now comes the trick |I], |, ^ |IU] of negative dimensions: Expand the integral ( |2.1| ) 
in Taylor series, 

^ _ ^ (-l)*+J+'+™a^/?jyc<j'" 

~ . n 7!7!/!7?7! 



X / I d'^qd^'r {q'y {q - pf \r^)' {r-q + kf"\ (2.5) 



and let us define 



Jndim = I d'^q d'^r (qj {q - pf ' (r^)' (r - g + kf . (2.6) 



We already note that this would be exactly the Feynman integral needed to be eval- 
uated for the diagram in Fig.l if it were not for the positive exponents 7, / and m. 
Comparing now (|2.4|) and (|2.5|) we conclude that, 



n^g{i,j,l,m) ~ (V)"^+"^r(l - m - n2 - \D) 

~ 1^^^ 77i!772!773!7l4!775!776! 777! ^ 



where 



(2.7) 

(7(7, J, /, 777) = r(i + 7)r(i + j)r(i + /)r(i + 777). 

We can immediately see that there is a seven index summation^ and five equations 
linking them (remembering that one comes from the constraint imposed by the multino- 
mial expansion). Therefore, altogether the system can be solved in twenty-one different 
ways with two remaining series^. Nine of them are trivial solutions, which present no 
interest at all. The remaining twelve we must solve one by one. In principle we do not 
know whether these are equivalent or not 0, [Tl|, |12| . 

A little bit of algebraic rearrangement yields, 

Jndim = {-ti)"" {p'T g{i. 3.l.rn)T{l ~ a -]^D) 

^ "ni+n2+n3+n4,iOni+n2+n5+n6,j"rii+n3+n5+n7,iOn2+n4+n6+n7,m /2 g\ 

l^'r'o} 77i!7l2!r73!774!775!776!n7! 



^i.e., a "heptaple" series 
^the remnant double series 
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where we have defined a = i + j + l + m + D. 

It is a simple matter to write down a computer program that solves the system exactly 
in all its twenty-one different ways. One of the solutions, namely, that with remaining 
sum indices 72,2 and uq is. 



Si 



-^)^(p^)"r(i + 2)r(i + j)r(i + /)r(i + m)r(i -a -Id) 

r(l + a)r(l - J - m - lD)T{l - I - lD)T{l + j - a) 

1 i-(r\n2){lD + l\n2 + ne){a-j\ne) 



(2.9) 



E 



where 



r(l + / + m + ID) n2W{l -j-m- \D\n2 + n^) 

T{a + h) 



[a\h) = {a)h 



Via) 



is the Pochhammer symbol ||13|, |15[ and we use one of its properties, i.e., 

(-1)' 



{a\ — A;) 



(l-a|A;) 



(2.10) 



within the double series. Note that the sum above can be rewritten in terms of hy- 
pergeometric functions 2F1 if we use another property of the Pochhammer symbol, i.e., 
{a\b + c) = (a + b\c){a\b). Now, a hypergeometric function 2-F1 with unit argument can, 
within certain constraints in its arguments, be summed |P^, |T1|, IT3|, 



2Fi(a,6;c|l) = 
Rearranging first the ng sum we have. 



T{c)T{c-a-b) 
T{c-a)T{c-b)' 



(2.11) 



(-7r)^(p^)-r(i + z)r(i + j)r(i + /)r(i + m)r(i -a-^D) ~ 

r(l + a)r(l -J-m - iZ})r(l - l - ID)T{1 +j-a) 

{-a\n2){\D + l\n2) ^ (a - 3\n^){\D + / + ^2^6) 
X — —. r 2^ 



n2!(l - j- m- \D\n2) 7161(1 - j - m-\D + n2\n&) 
1 

X— — -, 2.12 

r(l + / + m + iD)' ' 

where the second series is by definition the gaussian hypergeometric function, 2-^1. Using 
( p.ll|) we can sum it and then sum also the series in n2, to get. 



Si 



'(p')"r(l + i)T{l + ])T{1 + /)r(l + m)r(l -a-\D) 



r(i + a)r(i 
r(i-/- 



m 



\D)T{1 + J - a)r(l + l + m+\D) 
D) 



T{l + i-a)T{l 



m 



D) 



(2.13) 



5 



Now, grouping the gamma functions in convenient Pochhammer symbols and using 
( p.lO| ) we analytic continue to negative values of the exponents I, and m and we go 
back to positive D to get, 



S^^ = n''{pY{-^\a){-J\a)i-l\-m-^D)i-m\a-^-J-^D) 

x(a + ^D\ -2a- \d){D + l + m\-l- ]^D). (2.14) 

This is the general result, in Euclidean space, for the Feynman graph of Fig.l. 
In the important particular case when i = j = I = m = —1, the result is 



cAC _ -n^'ipT-'T^D - ^)T\\D - l)r(2 - \D)V{A - D) 

~ r(D-2)r(|D-4) ' ^^-^^^ 



which is the well-known result in D- dimensions [|16|. 

We ask immediately: What is the result that the other solutions provide? The answer 
is as surprising as it could be: the result is the same. All the twelve non-trivial solutions, 
even distinct, give the same result. It is an amazing feature revealed by NDIM. 

Just to check on this, let us consider another solution, for instance, the one where n2 
and 71,5 are the remaining sum indices. 



^2 



-7r)^(p2)-r(i + z)r(i + j)r(i + or(i + m)r(i -a-^D) 



r(l + a)T{l -a-m- ^D)T{1 + j- a)T{l + / + m + ^D) 
1 



2 

oo 

E 





n2)i-j + crln^) 


n2!n5!(l — a — m — 


n2 - 





1 



(1 + z + m + 2-D|r25 - 712) 

As in the previous case we can sum both series. Using (|2.10 ) one can rewrite 
series and identify it 2F1 function. 



S2 = {-n)''{p'rP2{i,j,l,m;D) ^ 



-lY-(-a\n2) 



n2 = 



E 



n2!(l -a-m- lD\n2) 
(cr - j\n5){a + m + - 7121^5) 



(2.16) 

the 77-5 

(2.17) 



where 



P2{i,j,l,m;D) 



1 + z + 777, + |D| - 772) n5!(l + i + m + - n2\n5) 

r(l + t)T{l + j)T{l + 1)T{1 + m)r(l -a- ID) 



r(l + a)r(l -a-m- ^D)T{1 + j - a) 
1 



X : 



r(i + / + 772 + iD)r(i + Z + m + |D) ' 

Summing the 7^5 series with formula ( p.ll| ) and using again ( |2.10| ) to rewrite the n2 
series, we get 
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= (-vr) (p) r(l + . - a)r(l - / - i^) 

^ {-(r\n2){\D + l\n2) 
X y — TT — — ' 2.18 



This series is by definition a summable 2-^1 function; witli tlie lielp of eq. (|2.11 
get tlie expression (|2.13|) whicli leads to the correct result ( p.l5|) . 



we 



The reader can prove, following the same procedure, that all the twelve solutions 
provide the correct result. The question that arises is: Why is this so? We have no 
answer to this puzzle at the moment and can only conjecture that maybe if the remaining 
series has unity argument and is summable then the result will be degenerate. Of course, 
further research is necessary in order to prove or disprove this conjecture. 



3 Conclusion. 



Our two-loop "lab testing" for NDIM approach to calculate Feynman integrals has re- 
vealed some very interesting features of the method. The methodology is quite sim- 
ple: solving gaussian integrals and systems of linear algebraic equations. NDIM yielded 
twelve non-trivial solutions which give the same result, eq.( |2.14| ), for the general case, 
i.e., D-dimensions and arbitrary exponents of propagators. This work encourages us to 
tackle a more difficult task: the calculation of massive four point one-loop integrals]^, [ill 



and off-shell two-loop Feynman graphs [12|. Work in this line is in progress, and we 



have already obtained some more encouraging results. For the latter, for example, a 
new surprise with NDIM yielding twenty-four distinct (and new!) results, some of them 
in terms of Appel's|jT^ hypergeometric functions F4 which are simpler than the usual 
dilogarithms pTSfl . These new results will be the subject addressed in our shortly forth- 
coming paper. 
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A The Remaining Solutions. 

For the sake of completeness, in this appendix we list the remaining ten non-trivial so- 
lutions since in the article proper we have explicitly shown only two of them, and that 
these two give the correct result. All of these can be summed and analytically continued 
to positive D with the same ideas we used in section 2. 
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r{l + a + l + lD)r{l-l-lD)r{l + i-a)r{l+j -a) ^ 
^ 1 (-J + (y\n6){\D + /|n4 + n6)(-i + a\ni) 



(-7r)^(p^)-r(i + ^)r(i + j)r(i + /)r(i + m)r(i - a - ip) 
r(i + j + i + iD)r(i + i + m+|D)r(i + j-a)r(i + / + m + iD) ^ ■ ' 

1 ^ {-lY^+''^i-i + a\n^){-j + a\n^) 



X 



E 



r(l + i - (j) „^^=o ?^4!?T'5!(l +i + ^ + \D\n4^ - nr,){l + i + m + iD|n5 - n^Y 

(-7r)^(p2)<^r(i + z)r(i + j)v{i + or(i + m)r(i - a - 



r(l + i + Z + |£>)r(l + J + m + \D)T{1 +j- a)r{l + l + m+^D) 

1 ~ (_l)n3+n6(_i + ^|„3)(__^- + ^|„g) 



(A.3) 



X 



E 



r{l + i- a) ^^^^0 n3!n6!(l + i + / + |L>|n6 - n3)(l + i + m + iL'Ina - r^e) ' 
(-vr)^(/)"r(l + ^)r(l + j)T{l + 1)T{1 + m)r(l -a -ID) 



r(l + a + m + iD)r(l - m - iD)r(l + i - a)T{l + j - a) 

^ 1 y. {-I + ^|^3)(-i + ^|n5)(|l^ + m\ns + ng) 



(A.4) 



-TT 



'^(p2)-r(i + ^)^(l + j)r(i + /)r(i + m)r(i - ^ - ^d) 



" r(l + (7)r(l - z - m - |i^)r(l + i - a)T{l -i-\D) ^^'^"^ 

1 ^ (-a|n2)(-z + (j\ni){\D + l\n2 + n^) 

^ T{l + l + m+\D) nM{l-i-m-\D\n2 + ni) ' 



(-7r)^(j9^)"r(i + i)v{i + j)r(i + /)r(i + m)r(i - a - |d) 

r(l + a)V{l + J + m + |i5)r(l + i - a)V{l - m - \D - a)T{l + l + m + lD) 

(-l)"^+"^(-(7|n2)H + a|n3) 



X 



n2^^=o ^2!?^3!(1 +j + m + lD\n3 - n2){l - m - \D - a\n2 - n^Y 

(-^)^(p^)"r(i + i)T{i + j)r(i + /)r(i + m)r(i - - 

r(l + a)T{l + 1 + 1 + \D)T{1 + j- a)r(l - I - \D - a)T{l + l + m + \D) 

^ {-ir-+^^>{-a\n,){-3 + a\n,) 

nii=o m'^^elll + i + / + |/^|n6 - ni)(l - / - - ajm - ng) ' 
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X 



-j^ oo 

r(l + / + m + iD) „^5=o 









- m 


ril + Tlr,) 


ni!n5!(l - j 


- / - ^D\ni + ns) 



s 
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(-7r)^(p^)"r(i + z)r(i + j)r(i + /)r(i + m)r(i - a - |d) 

r(l + a)r(l + i + / + ID)T{1 + i - a)T{l - I - \D - a)V{l + I + m + \D) 

^ g (-l)"^+"^(-a|ni)(-^ + a|n4) ^ ^^^^ 

ni^=o '^i''^4!(l + i + / + |£'|n4 - ni){l - l- \D - a\ni - n^Y 



^ (-7r)^(p^)-r(i + i)v{i + j)r(i + /)r(i + m)r(i - a - ip) 

r(i + a)r(i-z-/-iD)r(i + z-a)r(i-m-|D) ^ ■ ^ 

oo 

""ni + Z + m + iD) „^5=o 
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